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1 Introduction 



Recent progress in understanding of lower dimensional string theories is par- 
tially due to the advantage of discrete methods. It happened that the matrix 
quantum mechanics, despite its simplicity and solvability describes nonper- 
turbative aspects of string theories at d=l target space, i.e. c = 1 matter 
coupled to the two dimensional gravity. 

The standard combinatorial methods of matrix models, as well as the 
continuous approaches have appeared so far to be ineffective for higher di- 
mensional target space. 

On the other hand, it has been known for a long time that certain matrix 
field theories are completely integrable in 2d for an arbitrary size of the 
matrix field. One of the most representative integrable matix field theories 
is the principal chiral field (PCF) which describes a free field on a principal 
, say SU{N), manifold: 

S = ^Jd'xtT [d,g^d,g] (1) 

where g is an NxN unitary matrix. Its large N solution has been anticipated 
for a long time to follow from its finite solution |2|, ^, ^. 

Apart from being a model of Goldstone bosons whose strong interaction 
is entirely determined by the geometry of the manifold, this model has a long 
history of exploring its analogy to QCD and contour geometry of gauge 
fields |]. 

The large N limit of the model is of particular interest. It is conceiv- 
able that it describes a string theory in two physical dimensions due to the 
analogy between planar Feynman graphs and the world sheets of a string. 
Of course one should not take this analogy literally: in asymptotically free 
theory neither the coupling constant Aq nor a renormalized coupling is a cos- 
mological constant of a string [0: due to renormalons the contribution of 
even planar graphs grows factorially with the order. However, some signs of 
the stringy behaviour will be seen in the nonperturbative regime. 

In this paper we present the exact and explicit large N solution for the 
chiral field in two dimensions. 

The exact solution of the PCF (both ^-matrix and Bethe Ansatz equa- 
tions) was found in Ref. |^. It turns out that the spectrum of the, say. 
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SU (N) model contains — 1 kinds of massive particles. They form multi- 
plets belonging to the diagonal of the SU{N) ® SU{N) associated with all 
fundamental representations of the SU{N) algebra, namely, the vector rep- 
resentation and all antisymmetric tensors according to the Dynkin diagram. 
The spectrum of masses is 

= ^ • /^N (2) 

where / = 1,...,A^— lis the rank of a fundamental representation and m = mi 
is the mass of the vector particle. In the two-loop approximation it is 

1 /— 4vr , , 

m = A—=exp{— — ) (3j 
V Ao Ao 

where A is a cutoff. All particles are bound states of the vector particles. 
At large N we must distinguish two physically different situations: 
N ^ (X) but m = mi = fixed. This means that mi = Imi, so that the 
Z-th particle is not a bound state any more. It is decomposed into / vector 
particles. This suggests that the interaction vanishes in this limit and we 
end up with a free massive field. This limit is not of a particular interest. 

Below we consider another and the only physically interesting limit: 
oo but the heaviest mass mN/2 = fJ, of the largest antisymmetric tensor re- 
mains fixed. In this case the masses fuse so that the mass spectrum becomes 
continuous. The label running along the Dynkin diagram becomes a contin- 
uous parameter. The energy of the massive excitation with a momentum p 



will be therefore yp^ -|- mf ^ ^p^ -f where we introduced q = fJ^j^l- We 
observe that an extra dimension emerges from the matrix structure of the 
field. This means that at = oo particles do not form a discrete spectrum, 
so the theory ceases to be a theory of point like particles. 

We shall find the energy of the ground state as a function of "the Noether's" 
field by adding a term tr^HiQL + HrQr)/2 to the hamiltonian of the the- 
ory corresponding to the lagrangian (1). Here Ql = I d^xgd^g'^ and Qr = 
J d^xg~^dog are the Noether's left and right charges and -ffi?(L) = diag{hi, h2 — 
hi, Hn-i — — ^iv-i) is an element of the Cartan subalgebra. It 

amounts to introducing in eq.(|ID the covariant derivative 

Df^g = d,g - td^oiHtg + gHR)/2 (4) 
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instead of the usual derivative df^g. In what follows, we shall consider only 
the case Hl = Hr = H. 

Parameters hi play the role of chemical potentials for elementary particles 
of the model, so that the energy £{h) is the energy of the ground state 
with a symmetry of the Young tableau [l-^-A^i^ 2-^i--^^ ...N^'^-^] where Mi = 
—d/dhiS{h). The field H introduces an energy scale into the theory and 
gives a valuable physical information. 

The large N solution is explicit. To ease the references we state it now. 

We use a special direction of the field which excites all types of the par- 
ticles on equal footing: 

^ sin(f /) 

hi = h . (5) 

sm(^) 

We show that the energy of the ground state is expressed in terms of modified 
Bessel functions: 

f{h) ^ ^^{S{h) - S{0)) = -^BH,{B)K,{B) (6) 

where the parameter B is defined through 

777 

- = BK,{B) (7) 

The distribution of rapidities of physical particles will obey the simple 
semi-circle law with the support B. The parameter B defines the value of 
rapidity corresponding to the Fermi momentum of the fused particles. We 
shall see, that B gives the most natural definition of the renormalized ( 
running) coupling constant : 

m = ^ (8) 

The reader may find some similarity between B and the Fermi level of eigen- 
values of the matrix field in the c=l string theory ||^. 
The paper is organised as follows: 

In section 2 we review the derivation of the spectral Bethe ansatz integral 
equations for rapidities. 

In section 3 we solve the spectral equations in the large N limit. 

In section 4 we obtain a singular behaviour on the treshold ~ m of the 
spectrum. 
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In section 5 we show that our solution agrees with perturbative two loops 
calculations at large h/m and find exact value of the mass ( namely the ratio 

In section 6 we calculate all terms of perturbation theory in running 
coupling constant A and show their factorial growth. 



2 /S-matrix and Bethe-ansatz equations for 
any 

Perhaps the most economical way to obtain the Bethe-ansatz equations for 
the chiral field is the factorized bootstrap method 0, rather than direct 
diagonalization of the hamiltonian of the model The point is that the 

^-matrix can be easily found on the basis of some heuristic hypothesis. Let 
us give a sketch of this approach. 

1)S -matrix. The chiral field is renormalizable and asymptotically free 
05 mi- is invariant under the left-hand and right-hand group transfor- 
mations G ^ G and the action (|l|) is defined only by the Lie algebra of G, 
i.e. does not depend on the representation of the group G. Therefore it is 
natural to assume that the elementary particles are massive and belong to 
some fundamental representations of the diagonal of the G 0G, whereas the 



antiparticles form conjugated representations. The model is integrable Wl 
therefore the scattering is factorized. Under these assumptions the minimal 
S'-matrix (factorized scattering matrix with a minimal set of singularities) can 
be determined unambiguously. The complete proof of the minimal S'-matrix 
being the scattering matrix of the chiral field requires a more sophisticated 
technique (see e.g. |jl], 0]). 

It turns out that once we assume that there is a particle in some, say, l-th, 
fundamental representation, the factorized bootstrap tells us that there are 
particles in all — 1 fundamental representations. In fact, they are bound 
states of an arbitrary chosen representation. Therefore it is convenient to 
start from the vector particle. The factorized SU{N) ® SU{N) scattering 
matrix for vector particles is the tensor product of the SU{N) factorized 
vector S'-matrices S = X{6)S{6) S{6). Here ^ is a rapidity of a massive 
relativistic particle (p° = mcosh^,p^ = msinh6') and X{6) is the CDD- 
ambiguity factor which cannot be determined by the factorization, unitary 
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and crossing symmetry conditions. The SU{N) unitary, crossing invariant, 
factorized S'-matrix of vector particles is well known |12|. It is 



where is the projection operator onto symmetric (antisymmetric) states. 

The amplitude in the symmetric channel u{6) and the amplitude in the 
cross channel (particle- antiparticle scattering) t{9) = xi2-\iN-eJ{2i-K) ^(^^ ~ ^) 
obey the unitarity conditions t{6)t{—6) = u{6)u{—6) = 1. The minimal 
solution of these equations is 

u{e) = (10) 
^ ^ rfi + —)T — - —) 

Finaly the CDD-factor is chosen to cancel all double zeros and double poles 
on the physical sheet < ImO < vr: 



?7r ^ 



This is the S'-matrix of the vector particles. It has a pole on the physical 
sheet at 9b = 2TTi/N in the antisymmetric channel. It corresponds to the 
first bound state ( the second rank antisymmetric tensor) with a mass m2 = 
m sin(27r/A^)/ sin(7r/A^). The S'-matrix of these particles can be also found 
by tensoring the vector S'-matrix (the fusion procedure). It also has a pole in 
the antisymmetric channel, and so on. In this way the whole mass spectrum 
(0) can be generated. 

2)Bethe-Ansatz Equations. The thermodynamical properties of the model 
can be obtained by imposing the periodic boundary conditions. For an inte- 
grable problem they imply the balance of two particle scattering phases and 
a phase of a free motion between collisions. For the i-th particle with the 
momentum msinh^j, the periodic boundary conditions lead to the problem 
of the diagonalization of a product of two-particle S-matrices. Say, for a state 
with A/" vector particles in the box L we have 

exp(imLsinh^a) = Yl ^ap{6a — 6p) (12) 
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The eigenvalues of the operator in the r.h.s of the eq.(|l2]) can be found by 
the projection onto the symmetric subspace 

M 

exp(imLsinh6'Q) = expiicpiO ^ — 6 p)) (13) 

where exp{i(j){6)) = v?{9)X{9) . To obtain the Bethe-Ansatz equation for the 
state which contains all kinds of particles, say, for the state with the Young 
tableau \^\-^--'^'^ ^2^^~^^ , ...N^'^-^] one has to consider complex rapidities of 
the bound states -"strings" 6 6^ + 2rT:i / N , where 6^ is a rapidity of the /-th 
particle and r is an integer running between —1/2 and 1/2. Substituting this 
into eq. ([T3|) and multiplying equations over r we shall obtain the equations 
for the rapidities of the state which contains Mi particles of the kind I: 

N-l Mn 

exp(zLm;sinh^«) = n 11 exp{ict)in{e^^ - Of^)) (14) 

n=l a=l,/3^a 

where 

Me)= J2 (f)ie + 2rTT/N + 2r'7T/N) (15) 

-l/2<r<l/2,~n/2<r'<n/2 

After tedious calculations 01 we obtain 



d9 
where 



oo 



dMd) ^ _2 / du[Rin{ij) - Sin] cosuO (16) 



sinh inuj(l — -^max(L n))) sinh (nuj-^min(L n) 

Riniuj) = 2 ^- ^ —-^ (17) 

smh TTU 



Taking logarithm of the both sides of the Eqs.([l^) we obtain the Bethe - 
Ansatz equations 

n 1 N-l N„ 

m.sinh^« = ^j« + |^ ^ 0.n(e«-<^) (18) 

where integers J are the quantum numbers of the states. For the generaliza- 
tion of these equations to an arbitrary Bethe states see 0, Q . The energy of 
this state is obviously 

■1 AT-l Ml 

^ = y E^'E^osh^W (19) 

^ (=1 a=l 
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3) Spectral Equations. The next step is to find rapidities to minimize 
the energy ([T9|) in the thermodynamic hmit Mi/L = ni, while L oo. We 
assume that the minimum of the energy corresponds to a dense smooth set 
of 6''s, so one can describe them by the distribution function of rapidities of 
particles pi{9) and the distribution of holes pi{0). They are related by the 
equation 

cosh^ = pi{e) + f de'Rinie - e')p^{e') (20) 

where 

1 /■°° 

Ri^(0) = - dujRin{uj) COS ujO (21) 

TT JO 

Let us now turn to the energy of the state at the given fields hi 

8 = mirinXE - ^i^i) = 51 / de[micosh6 - hi]pi (22) 

I •' 

Consider a small variation of piiO) and pi{6) and introduce the (pseudo) energy 
ef{6) > of a hole and the (pseudo) energy eY{6) < of a particle [Q, such 
as 5£ = ! d6[e^ {6)5pi{6) — {6)5pi{6)]. According to this definition at 
the ground state ei{e) = if pi{0) ^ and eY{e) = if pi{e) ^ 0, i.e. 
e^{ej) and pi{pi) have nonoverlapping supports. Comparing with (|2^) and 
using (^) we find the spectral equations of the model 

hi - mi cosh^ = ei{e) + Rin{e - e')e+{e'), (23) 

^ = -^E/^/^(^Hcosh^, (24) 

where e^O) > 0, er{9) < 0, e+{0)er{9) = 0. 

4) Diagonalization of the Spectral Equations. To prepare the spectral equa- 
tion (|23D for the large N limit, we diagonalize the kernel matrix Rin. This is 
easy to do since it reflects the structure of the Dynkin diagram @] . Indeed, 
its inverse is 

Rr„\uj) = Sinh-^ ^{^In cosh ^ - l/2{6l,n+l + 5i+l,n)) (25) 
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Let us introduce orthonormal eigenfunctions of the (A^ — l)x(A^ — 1) Cartan 
matrix Cjh = 25 jk — — 5j+i,fc: 



= y27iVsin^, p=l,2,...,iV-l 
2^k=i Xk Xk — '■'pp 



(26) 



They are also the eigenvectors of Ri 

RA-^) = E xfx'^'R^'Hu:) (27) 

where 



,(p)^ip) 

p=l 



p(p)/ N sinh^ _2N ^ \u\ 

^""^ - cosh ^- cos f ~ V ^hoo u^ + ip + rNy ^'^^^ 

Then the spectral equations can be diagonalized with respect to the particle 
indices: 



{6) + / R^p\e - e')ef{e')de' = h^p^ - 5p,iMcoshe, (29) 



^iP), 



£ = ^ ' / d^cosh^A (^) ^ - ^ . / d9 cosh 9e^!:\9) 

27rsin(7r/iV) J + ^ ^ ir^ J + ^ ^ 

(30) 

where 
and 

Af-l N~l 

^ E X?^4^^; h^'^ ^ E x?^/^. (32) 
fc=i fc=i 



It is important to note that this transformation is valid and eq.(p9D holds, 
providing that linear combination of particle's (hole's) (pseudo)energies re- 
main positive (negative): > {^'f'l) < 0). 
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3 Large solution 



One has to be careful using the last definition of the Fourrier transform of 
e±\uj) = /f^ d6e±\6) cos{uj6): different e^{6) entering the sum in (|32|) have 
different supports {—Bk,Bk). It has to be taken into account when solving 
the integral equations. 

However, if we take hk in the form (||), so that 



h^"^' = h5r, » Af^oo r- hSr, 1 (33) 

'^'V2sin(7r/Ar) V27i ^ ' 



it turns out that the ansatz 



ef\uj) = e±{uj)6p^i— - ^n^oo -^e±{uj)6p^i (34) 



^A^sin(7r/iV) 

satisfies all the equations (^) for p = 2,3, N—1 and the equation with p = 
1 (the Perron- Frobenius mode) gives an integral equation for the definition 
of e±{uj). Since e+{9) is strictly positive it can be viewed as a newly defined 
density, with all appropriate analytical properties. 

It implies that we look for a solution of (^) with equal supports Bi = 
B2 = ... = -Bat-i = B. In what follows we consider eq.(^) only inside the 
interval {—B, B) where e-{6) = 0, e^{9) = e{9) > 0. It is the only function 
which contributes to the ground state energy (pOf). It obeys eq.(29) with 
p = 1. 

Further simplifications occur in the large limit. The kernels R^^^ look 

as: 

2N \uj\ 



Rip)(u;) = iV^oo (35) 

Finally, for the choice ( p3D of the field we obtain the integral equation: 

/OO (^i^ I ijj I 

— cos(gLj) ' ' e(uj) = h~mcoshe, \0\ < B (36) 
^00 27r to'^ + 1 



where ^ 

e(uj) = I deeie) cosiOuj) 

J-B 

Note that, in the large limit R^p\uj) vanishes at large u, whereas at finite 
it approaches 1 (see eq.(28)). This implies that e{d) is no longer differentiable 
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at 6 = ±5, but has a cusp. As a result the physics on the threshold h ^ m 
will be changed drastically. 

Equation (BBI) can be solved exactly. Note that if we act by the operator 



(— ^ + 1) on both sides of it we obtain a simple integral equation with the 
singular kernel: 

n^J_si9-9'y~ " ^^^^ 
Its solution is the famous semi-circle law of Dyson: 



e{e) = hVB^ - 02 (38) 

Its Fourier transform can be expressed through the Bessel function 

e{uj) =nhBJi{Buj)/u (39) 

Plugging it into (^) and doing the exact integration over u (see [|l3l) we 
obtain the relation (]^) between the m/h and B. 

Now we can calculate the free energy as a function of B using the relation 
( |30D and our solution (|38|): 

f(h) = ^ (£(h) - 8(0)) = — \;mh r dO cosh O^B^ - 9^ (40) 

The integral (^O]) together with the relation (0) gives the result (||) for the 
energy . 

Let us note that in the large limit any virtual and real processes involve 
all particles, since minimal energies are greater then a minimal separation 
between masses. A reasonable external field (^ excites all of them on equal 
footing and leads to the collective effects. That is why, even though the 
S-matrix of two, say, vector particles tends to unity at A^ ^ 00, the ground 
state energy is not that of a free field theory. 



4 Singular behaviour on threshold 

Here we will show that the theory exhibits qualitatively new features on the 
threshold h —>■ m which corresponds to 5 — > 0. 
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At small B asymptotics of the Bessel functions Ii{B) 1/'k{B/2 + 
5716 + ...) and Ki{B) 1/5 + f ln(5/2) + .... Then, from (|) we obtain: 

B^ ~ 4—^, A ^ (41) 
|lnA| ^ ' 

where we introduced A = h/m — 1. This gives a singular behaviour on the 
threshold: 

^ -(m/27r)2^^, A^O (42) 

It differs drastically from the threshold behaviour for a finite N theory of 
massive particles, where we would have 3/2 law (see e.g. 0, 



Mh) ~ -m\A)'/' (43) 

As we already mentioned, the reason for the singular behaviour is the 
emergency of an extra dimension in the large limit - the masses of physical 
particles are distanced by each other by the quantity of the order 1/A^, which 
is less than any energy scale left in the system. Therefore any external field 
excites a bundle of particles - a new extended object, which is characterized 
by an extra "momentum" q in addition to the usual momentum p. Let us 
note that a similar behaviour has been found in the quasiclassical limit of 
the Sine-Gordon model |T^, [0. 



It is interesting to compare the results (|41|) and ( ^21 ) with the expression 
for the ground state energy of the c = 1 matrix model ||^. The mechanism 
by which the inverse logarithmic behaviour with respect to the cosmological 
constant occurs also requires a parametrization through the fermi level of the 
corresponding fermions (whose coordinates are the eigenvalues of a hermitean 
matrix field). The Fermi level plays the role of a "hidden" parameter of the 
problem and the eigenvalues give rise to an extra (Liouville) degree of freedom 
of the theory [|17 . 



5 Perturbative regime and the ratio m/ K^g 

The field h fixes the scale of energies (with respect to the mass gap m). 
Large h/m regime is the subject of the perturbative theory. From eq. (|^) 
we conclude that it corresponds to i? — >^ cxd. It follows from the large B 
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asymptotics of McDonald function Ki{B) = ^f^e'^yl + ^ + 0(1/5^ 
that 

h [2 3 



m V TT ^/B ^ SB 
Solving for B we obtain: 



1- — + 0(1/5^) (44) 



B = \n- + -ln\n- + -ln- + 0(^) (45) 
m 2 m 2 2 Mn>^ ^ ^ 



Using the large B asymptotics Ii{B) = ^^e^ (l — -^ + 0(1/ B"^)^ we finally 
find from flBD and fHSl): 



IQnfih) = -h'B+0& = -h'(\n- + l\n\n- + l-\n^ + 0{-^)) (46) 

B ^ m 2 m 2 2 In— ^ 



This result reproduces correctly one- and two-loop terms of the perturba- 
tion theory as well as the universal non-perturbative constant | In ^ (first 
calculation of a similar constant was given in |19] for the Kondo problem). 

Let us compare it with the perturbative result following directly from 
(|l]). In our normalizations we have 

h"^ ^ 2h^ 1 h'^ 

(47) 

where 

qk = ^'~j'^''\ k = l,2,...N; ho = = (48) 
and the renormalized coupling in the minimal subtraction {MS) scheme is: 

(/^) = ^ ( In + ^ In In + 0(1/ In -^)) (49) 

4vr^ 2Ams 2 Ams Aa75 ^ 

Or, 

,|-A(/,).^(A) = -^P-J-,Ph-0(A') (50) 
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Using an approximate solution of the BA equations at finite for a particular 
choice of the field which excites only the vector particles with the lowest mass 
mi, the authors of found the following ratio: 

m r — — sin(7r/A^) 



/ — Sm TT/JV / — , , 

A \We .1/^' ^N^oo (51 

To compare it with our result ( ^Bj ) we have to calculate all the sums in (^) 
for our choice of fields (|]) : 

cos jf{k — i) 

qk = ^ ^AT-^ooCOSX, < X < TT (52) 

cos 

The calculation of corresponding sums in the large A^ limit gives: 

^Ell-^N r^cos^x = 1/2 (53) 

W T.k>Mk - gi)^[ln \qk - gil - 1/2] 

1/2 Jo" f Jo" f (cosx - cosy)2[ln | cosx - cosy| - \] = {I - 21n2)/4 

(54) 

Now matching our result (^61) with the perturbation theory (^T]) we obtain 
the same ratio (|5lD . So, two different choices of the field (that of |20| and of 
our's) give the same mass gap in the MS scheme. 

Let us note in the conclusion to this section that we can continue compar- 
ing further terms of the expansion ( ^Bj ) of our exact ground state energy in 
the inverse logarithms with the planar part of the Feynman graph expansion 
( ^7D for the ground state energy using relation (^Tj). They must coincide in 



all orders. 



6 Weak coupling expansion to any order and 
exact beta-function 



The \{h) defined eq.(ES|) is equal to An / B in the two-loop approximation 



(^5]). Since all next corrections are nonuniversal we may take An/B as the 
most natural definition of the renormalized coupling: 

X{h) = Att/B (55) 
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Having in hands the exphcit expression for the free energy we can 
find all orders of the renormalized perturbation theory. Let us notice that 
the function: 

y{B) = BK^{B)h{B) (56) 
satisfies the 3-d order differential equation: 

l/'"-(j^ + %'+|^^ = (57) 

If we look for the series expansion of y{B) we find from (0) the recurrence 
relations on the coefficients, which can be solved. Finally, by taking (^) into 
account we obtain: 

/(A)/A' = -i(l-|c.„(A)-) (58) 



where 



2n + l((2n)!)' ^ ^,n.,^ 



Note that the expansion goes only in even powers. 

As we see, in spite of the fact that every coefficient represents a sum 
over renormalized planar graphs, it grows factorially with the order. Most 
probably this happens because of the renormalons (some subsequence of log- 
arithmically divergent graphs) giving the main factorial contribution in each 
order noticed long time ago by 'tHooft 0]. This means that we have an 
exponential number of graphs in each order but some of them give (2n)! 
contribution after the momenta integration. More than that: the series is 
a non-signchanging one and thus non-Borel summable. Nevertheless, the 
free energy perfectly exists for any finite A. These phenomena seem to be 
imminent for any asymptotically free field theory. 

One can still give a prescription how to sum up the series ( p^D and restore 
the result (1). Notice that: 

oo 

9it) = 1 - E 7^^'" = F{-l/2, 3/2, 1, tV4) (60) 

n=l 

where F is the hypergeometric function. It has the cut starting from = 4. 
We can continue g{t) analytically to the whole complex plane and make the 
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inverse Borel transform by integrating from to oo with the exponential 
factor. The correct prescription restoring (Q) is (see (|55|)): 



16nf{h)/h'^ = -B^Re dte-'^g{t) (61) 

Jo 

As we see from here even the Borel nonsummable series in the asymp- 
totically free theories can be summed up by some special prescription. This 
prescription which uses a non-trivial procedure of analytical continuation 
might work well in other asymptotically free theories. Note that the inte- 
gral in ( |6TD taken along the real axes possesses also an exponentially small 
( for B oo, A — >■ 0) imaginary part equal to —K1{B)B'^/tt. The similar 
properties of Borel integrals in 4d gauge theories were noticed in 

With the definition (|^) of the running charge one can find from eq. (0) 
the exact beta-function: 



or 



where 



bo = 1, bi= 4, 62 = -8, 63 = 64, 64 = ■ 2^ 
65 = 13 ■ 210, &6 = -1073 • 2^, , 67 = 103 • 216, 



(64) 

6„ ~ -{-ir^8/inn){An/er (65) 



7 Discussion 

Our results could be interesting from several physical points of view: 

l.The model considered here is an example of an exactly solvable matrix 
model in 1+1 physical dimensions. A temptation would be to interprete it as 
a new string theory in 1+1 dimensional target space. Recall that a natural 
interpretation of big planar graphs could be given in terms of a sum of world 
sheets of a string (see and references therein). An important condition for 
it is that the typical graphs should be big and dense, to describe a smooth 
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surface. For this one usually tunes the coupling to its (finite) critical value, 
corresponding to the "explosion" of the size of graphs. However, in the 
asympthotically free theory the critical coupling is equal to zero. It is the 
consequence of the fact that the theory has exponential corrections to the 
perturbative expansion, leading to the factorial divergency of its coefficients. 
In planar theories it can happen because of the presence of renormalons: some 
small portion of all the graphs has a factorially big weight (with respect to the 
order). These graphs seem to be not very much suitable for the interpretation 
in terms of random surfaces. Even the expansion in terms of the renormalized 
coupling contains the same factorial contributions, as we see from our results. 
The same scenario should be available to the multicolour QCD, invalidating 
this naive relation between planar graphs and world sheets of a hypothetic 
QCD string. 

2. However, for both PCF and QCD a less trivial string scenario could be 
possible. It might exist a (nonperturbative) reformulation in terms of a new 
master field, with a new perturbation theory, already suitable for a string 
interpretation. The result ( ^2]) suspiciously reminds the similar result for 
the ID bosonic string (from the matrix quantum mechanics) [§, if we take 
{h — m)^ instead of the cosmological constant A — Xcrit- The similarity is 
even more striking if one remembers that for both models this asymptotics is 
related to the behaviour of the fermions on the top of the fermi-sea (|]1], 0), 
and the Fermi level is the most natural physical parameter in both cases. 
The semi-circle law for the distribution of rapidities (37) also suggests that 
at some value of field h we deal with an extended object. 

3. Another interesting feature of the model (also observed in the c = 1 
matrix model) is the emergency of an extra dimension following from the 
matrix structure of the theories. This dimension is related to the random 
walk along A^^i Dynkin diagram. The structure of this extra (third) dimen- 
sion is well seen from the fact that the kernel of the problem (see eqs. (|25|) 
or (pSD) looks as a propagator of a periodic motion in the space of rapidities 
versus Dynkin diagram. It reminds an extra (eigenvalue) dimension in the 
collective field formulation of the ID matrix model usually attributed to the 
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4. One of the immediate problems is to take into account the effects of 
finite N. The most interesting effects of finite A^'s correspond to the e~^^ 
corrections. To demonstrate it let us look at the expansion of the kernel (^) 
in the pole terms (see (|2^)). The main order for the large N corresponds to 
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the neglection of all poles but the first (on the finite distance from the origin 
on the imaginary axes of a;). Next terms will give the exponentially small 
contributions (as well as the 1/N corrections) if we treat them as singularities 
in the corresponding integral equation. 
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